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Abstract 

The behavior of highly charged short rod-like polyelectrolytes near oppositely charged planar 
surfaces is investigated by means of Monte Carlo simulations. A detailed microstructural study, 
including monomer and fluid charge distribution, and chain orientation, is provided. The influence 
of chain length, substrate's surface-charge-density and image forces is considered. Due to the lower 
chain-entropy (compared to flexible chains), our simulation data show that rod-like polyelectrolytes 
can, in general, better adsorb than flexible ones do. Nonetheless, at low substrate-dielectric- 
constant, it is found that repulsive image forces tend to significantly reduce this discrepancy. 
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I. INTRODUCTION 



Polyelectrolyte (PE - i.e., charged polymer) adsorption on charged surfaces is a versa- 
tile process having industrial as well as biological applications. In particular, the case of 
rod-like (stiff) PEs, which corresponds to the situation of short DNA fragments or other 
similar biomaterials, has some relevance for biological processes such as gene therapy^ or 
biotechnology.—^ 

From a theoretical viewpoint, the behavior of stiff PEs near an oppositely charged surface 
has been studied by various authors on a mean-field leve L^ 6 ' 7 ^ 9 ' 1 ^ 11 Menes et al.— considered 
the interaction between two infinitely long charged rods near a salty surface in the framework 
of the Debye-Huckel theory. Due to the low dimensionality of the system they reported an 
algebraic decay for the effective interaction that was confirmed by their Brownian dynamics 
simulations.— The more simple and fundamental situation of a single and infinitely long 
charged rod in the vicinity of a charged surface was investigated by several authors.^^ The 
problem of the so-called "counterion release" (i.e., "Manning decondensation") from a rigid 
PE approaching an oppositely charged was examined by Sens and Joanny 7 - and by Fleck and 
von Grunberg^ using Poisson-Boltzmann Theory. In a similar spirit, Menes et al.— found 
that the screening of the adsorbed charged rod's field, due to counterions and mobile surface 
charges, is highly sensitive to the degree of membrane's surface charging. The adsorption of 
short rod-like PEs was also considered by some researchers.—^ Recently, Cheng and Olvera 
de la Cruz^ investigated the adsorption/desorption transition including lateral correlations 
by assuming a regular flat lattice for the adsorbed charged rods. This latter assumption is 
only valid when the electrostatic rod-surface binding is strong enough. Closely related to 
our problem, Hoagland 10 analyzed the monomer concentration profile as well as the chain- 
orientation with respect to the charged substrate's surface for a single short rod-like PE. The 
notorious complication of image forces stemming from the dielectric discontinuity between 
the substrate and the solvent was also addressed by some authors . 7111 ' 12 More specifically, for 
low dielectric constant (i.e., repulsive image forces) and fixed surface-ions of the substrate: 
(i) Sens and Joanny 7 - showed that the condensed counterions are not always released as the 
stiff PE approaches the substrate and (ii) Netz^i showed an extra decrease of the charge 
dissociation of the stiff PE (i.e., charge regulation in salty solution nearby an uncharged 
polarized interface). 
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Since those studies ^ 10 ' 11 were realized in the framework of the Poisson-Boltzmann theory 
and for a single chain, the relevant phenomenon of charge reversal of the substrate's surface- 
charge by the adsorbed PEs can not be captured. 

In this paper, we propose to elucidate the microstructural behavior of (very) short rod-like 
PEs near an oppositely charged surface by using Monte Carlo (MC) computer simulations. 
The effect of image forces is also systematically investigated. To better understand the effect 
of chain-entropy, a comparison with the previous work of Messina^ concerning fully flexible 
PEs is carried out. Our article is organized as follows: The simulation model is detailed in 
Sec. HU Our results are presented in Sec. and concluding remarks are provided in Sec. 

hyi 

II. MODEL AND PARAMETERS 
A. Simulation model 

The model system under consideration is similar to that recently investigated for the 
adsorption of flexible chains . 13 ' 14 Within the framework of the primitive model we consider 
a PE solution near a charged hard wall with an implicit solvent of relative permittivity 
tsoiv ~ 80 (i.e., water at z > 0). The substrate located at z < is characterized by a relative 
permittivity e su b s which leads to a dielectric jump A e (when e so ;„ ^ e su b s ) at the interface 
(positioned at z — 0) defined as 

^ ^solv ^-subs /j\ 

(■solv ^-subs 

The negative bare surface- charge density of the substrate's interface is — o"oe, where e is the 
(positive) elementary charge and do > is the number of charges per unit area. The latter is 
always electrically compensated by its accompanying monovalent counterions of charge Z + e 
(i.e., monovalent cations with Z + = +1) of diameter a. Rod-like PE chains are made up of 
N m monovalent positively charged monomers (Z m = Z + = +1) of diameter a. The bond 
length I is also set to I = a so that the length L ro( i of a rod-like PE is L ro d = N m l = N m a. 
The counterions (monovalent anions: Z_ = —1) of the PEs are also explicitly taken into 
account with the same parameters, up to the charge-sign, as the substrate's counterions. 
Hence, all the constitutive microions are monovalent [Z = Z + = Z m = —Z_ = 1) and 
monosized with diameter a. All these particles are immersed in a rectangular L x L x r 
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box. Periodic boundary conditions are applied in the (x, y) directions, whereas hard walls 
are present at z = (location of the charged interface) and z = r (location of an uncharged 
wall). It is to say that we work in the framework of the cell model. 
The total energy of interaction of the system can be written as 



tot 



i.Kj 



-(intf) 
hs 1 



Zi) + U 



(intf) 



coul 



+ 



(2) 



(mic) 
hs I 



. TT (mic), 



where the first (single) sum stems from the interaction between a microion i [located at z = zi 
with i = (+, —,m)} and the charged interface, and the second (double) sum stems from the 
pair interaction between micro ions i and j with = |r^ — r 3 -|. All these contributions to 
U to t in Eq. (j2J are described in detail below. 

Excluded volume interactions are modeled via a hardcore potential defined as follows 
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oo, otherwise 
for the interface- micro ion one. 

The electrostatic energy of interaction between two microions i and j reads 
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where Ib = fie 2 /(47reoe so /„) is the Bjerrum length corresponding to the distance at which 
two protonic charges interact with 1/(3 = ksT, and A e is given by Eq. The first term in 
Eq. (j3J) corresponds to the direct Coulomb interaction between real microions, whereas the 
second term represents the interaction between the real microion i and the image of microion 
j. By symmetry, the latter also describes the interaction between the real microion j and 
the image of microion i yielding an implicit prefactor 1/2 in Eq.(jSJ). The electrostatic energy 
of interaction between a microion i and the (uniformly) charged interface reads 
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The second term in Eq.© stands for the self-image interaction, i.e., the interaction between 
the microion i and its own image. An appropriate and efficient modified Lekner sum was 
utilized to compute the electrostatic interactions with periodicity in two directions.— This 
latter technique was already successfully applied to the case of PE multilayering^ and 
polycation adsorption.— To link our simulation parameters to experimental units and room 
temperature (T = 298K) we choose a = 4.25 A leading to the Bjerrum length of water Ib = 
1.68a = 7.14 A. In order to investigate the effect of image forces we take a value of e su f, s = 2 
for the dielectric constant of the charged substrate (which is a typical value for silica or 
mica substrates^) and e so i v = 80 for that of the aqueous solvent yielding A t = ps 0.951. 
The case of identical dielectric constants e su b s = e S oiv (A £ = 0) corresponds to the situation 
where there are no image charges. 

All the simulation parameters are gathered in Table |U The set of simulated systems can 
be found in Table |H] The equilibrium properties of our model system were obtained by 
using standard canonical MC simulations following the Metropolis scheme . 17118 In detail: 

• Single particle (translational) moves were applied to the counterions (i.e., anions and 
cations) with an acceptance ratio of 50%. 

• As far as trial moves for the rod-like PEs are concerned and given the anisotropy of 



TABLE I: List of key parameters with some fixed values. 



Parameters 



T = 298K room temperature 

uqI? charge number of the substrate 

A e = or 0.951 dielectric discontinuity 

Z = 1 microion valence 

a = 4.25 A microion diameter 

Ib = 1.68a = 7.14 A Bjerrum length 

L = 25a (x, y)-box length 

t = 75a z-box length 

N ro d number of rod-like PEs 

N m number of monomers per rod-like chain 
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TABLE II: Simulated systems' parameters. The number of counterions (cations and anions) 
ensuring the overall electroneutrality of the system is not indicated. 



System N rod N m a L 2 

A 96 2 64 

B 48 4 64 

C 32 6 64 

D 24 8 64 

E 16 12 64 

F 24 8 32 

G 24 8 128 



these objects, random translational moves as well as rotational ones were performed at 
the same frequency.— Random rotational moves were achieved by choosing randomly 
new orientation-vectors of the rod-like particles. This method is (i) computationally 
not too demanding, (ii) leads to an efficient configurational space sampling and (iii) 
fulfills the rules of detailed-balance. The acceptance ratio was also set to 50%. 

The total length of a simulation run is set to 3 x 10 6 MC steps per particle. Typically, about 
10 5 MC steps were required for equilibration, and 2.5 x 10 6 MC steps were used to perform 
measurements. 



B. Measured quantities 

We briefly describe the different observables that are going to be measured. In order 
to study the PE adsorption, we compute the monomer density n(z) that is normalized as 
follows 

T-a/2 

n(z)L 2 dz = N rod N m . (7) 

a/2 

To further characterize the PE adsorption, we also compute the total number of accumulated 
monomers N(z) within a distance z from the charged interface that is given by 

N(z) = [ n(z')L 2 dz'. (8) 

J a/2 



It is useful to introduce the fraction of adsorbed monomers, N*(z), which is defined as 
follows 

The orientation of the rod-like PEs can be best monitored by the angle 9 formed between 
the z-axis and the PE-axis.-^ A convenient quantity is provided by its second order Legendre 
polynomial: 

P 2 [cos0(z)] = l - [3 cos 2 #(2) - 1] , (10) 
where z corresponds to the smallest wall-monomer distance for a given PE. Thereby 

S(z) = (P 2 [cos9(z)]) (11) 

takes the values — ~, 0, and +1 for PEs that are perpendicular, randomly oriented, and 
parallel to the z-axis, respectively. 

Another relevant quantity is the global net fluid charge a(z) which reads 



a(z) = / [n + {z') -n_{z')]dz', (12) 

J a/2 

where n + and n_ stand for the density of all the positive microions (i.e., monomers and 
substrate's counterions) and negative microions (i.e., PEs' counterions) , respectively. The 
corresponding reduced surface charge density o~*{z) is given by: 

a (z) = ——. (13) 
0q 

Thereby, cr*(z) corresponds, up to a prefactor a e, to the net fluid charge per unit area 
(omitting the surface charge density — o§e of the substrate) within a distance z from the 
charged wall. At the uncharged wall, electroneutrality imposes a*(z = t — a/2) = 1. By 
simple application of the Gauss' law, [a* (z) — 1] is directly proportional to the mean electric 
field at z. Therefore c*(z) can measure the screening strength of the substrate by the 
neighboring solute charged species. 

III. RESULTS AND DISCUSSION 

It is well known that the effects of image forces become especially relevant at low surface 
charge density of the interface. Furthermore, it is also clear that the self-image interaction 
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(repulsive for A e > 0, as is presently the case) is higher the higher the charge of the ions 
(polyions). Hence, we are going to study (i) the influence of chain length (Sec. IIII A|) and 
(ii) that of surface charge density (Sec. 1111 B|) . For the sake of consistency, we set the total 
number of monomers to N TO( iN m = 192 meaning that the monomer concentration is fixed 
leading to a PE volume fraction = | ^W) 3 ^ 2.14 x 10" 3 (see also Table EJ). 

A. Influence of chain length 

In this part, we consider the influence of chain length N m at fixed surface-charge-density 
parameter o§I? = 64. The latter would experimentally correspond to a moderate^ surface 
charge density with — o§e ~ —0.091 C/m 2 . The chain length is varied from N m = 2 up to 
N m = 12 (systems A — E, see Table ITT]) . We have ensured that, for the longest chains with 
N m = 12, finite size effects are not important since there L ro d = 12a which is significantly 
smaller than L = 25a or r = 75a. 

1. Monomer distribution 

The profiles of the monomer distribution n(z) can be found in Fig. Qand the correspond- 
ing microstructures are sketched in Fig. El When no image charges are present [A e = - 
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FIG. 1: Profiles of the monomer density n(z) for different chain length N m with o~qL 2 = 64 
(systems A - E). (a) A e = 0. (b) A e = 0.951. 
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FIG. 2: Typical equilibrium micro-structures of systems A — E. The little counterions are omitted 
for clarity. 
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Fig. IB^a)], the monomer- density profile n(z) exhibits a monotonic behavior for very short 
rigid chains (here N m = 2). For longer chains there exits a small monomer depletion near 
the charged wall for an intermediate regime of N m (here 4 < N m < 8). At high enough 
enough N m (here N m = 12) our simulation data reveal again a monotonic behavior of n(z). 
This interesting behavior is the results of two antagonistic driving forces, namely (i) chain- 
entropy and (ii) the electrostatic wall-monomer attraction. More precisely, the mechanisms 
responsible for this iV m - induced reentrant behavior at A e = are as follows: 

• For very short chains (here N m = 2) chain-entropy effects are negligible so that one 
gets a similar behavior to that of point-like counterions with the usual monotonic 
decaying n( 2:) -profile. 

• The chain-entropy loss (per chain) by adsorption should typically scale like lnN m 
whereas the electrostatic wall-chain attraction scales like N m explaining why at high 
enough N m a purely effective attractive wall-monomer interaction is recovered. 

Upon polarizing the charged interface [A £ = 0.951 - Fig. IHb)] the PE adsorption becomes 
weaker and the n(z)-profile more broadened due to the repulsive image-polyion interaction. 
For N m < 4, n(z) presents a maximum at z = z* pa 0.9a that can be seen as the thick- 
ness of the adsorbed PE layer. Interestingly, the monomer density at contact decreases 
with increasing N m . This is the result of a combined effect of (i) chain-entropy loss near 
the interface and (ii) the iV m -induced image-polyion repulsion. All those features are well 
illustrated on the microstructures of Fig. |21 

It is instructive to compare the above findings with those obtained for fully flexible chains. 
To do so, we use existing MC data for flexible chains from our previous worki^ where all 
the parameters, up to the chain flexibility, are identical with those presently employed for 
rod-like PEs. The comparison is provided in Fig. El At A e = [see Fig. Efa)], the n(z)- 
profiles for flexible and rigid PEs are quasi-identical for N m = 2 as it should be. For longer 
chains (N m = 8), we clearly see at A e = that the degree of adsorption as indicated by the 
value of n(z) near contact is considerably stronger for rigid chains. This feature is due to 
entropy and electrostatic effects. Indeed, in the bulk and at given degree of polymerization 
N m , the chain-entropy associated to rigid PEs is much lower than that associated to flexible 
chains, so that chain-entropy loss upon adsorption is reduced for rigid chains. Secondly and 
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FIG. 3: Comparison between flexible and rod-like (rigid) PEs for the profiles of the monomer 
density n(z) at different chain length N m . (a) A e = 0. (b) A e = 0.951. 

concomitantly, the wall-PE attraction is more efficient for rigid chains than for flexible chains 
because in the latter case the z- fluctuations of the charged monomers are more important. 

As far as the monomer density at true contact is concerned [i.e., n(z — > a/2)], it seems 
that, for flexible PEs, its value is nearly independent of N m as already reported in Ref.— . For 
a single ionic species of spherical shape, a variant of the contact theorem provides the exact 
relation: n(a/2) — n(r — a/2) = 2-ko\Ib yielding to n(a/2) pa 0.11a -3 , which is surprisingly 
in remarkable agreement with the value reported in Fig. [31(a) for flexible PE (and rigid ones 
for N m < 4). Nonetheless, Fig. Efa) shows, already with N m = 8, a strong deviation from 
the contact theorem (which in principle only holds for structureless spherical ions) for rigid 
PEs, as expected. 

The scenario becomes qualitatively different when A e = 0.951 [see Fig. Efb)], where the 
n(z)-profiles for flexible and rigid PEs become more similar. It is to say that the image- 
polyion repulsion tends to cancel chain-entropy effects. A closer look at Fig. Efb) reveals 
however that, at given N m , the degree of PE adsorption is systematically larger for rigid PEs 
than for flexible ones as expected. Those relevant findings can be summarized as follows: 

• Without dielectric discontinuity (A e = 0) rigid PE chains can much better adsorb 
than flexible ones at oppositely charged surfaces essentially because of their significant 
lower chain-entropy. 
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FIG. 4: Profiles of the fraction of adsorbed monomers N*(z) for different chain length N m with 
aoL 2 = 64 (systems A — E). (a) A e = 0. (b) A e = 0.951. The inset is merely a magnification near 
contact. 

• In the presence of polarization charges (A e = 0.951) the degree of PE adsorption 
becomes significantly less sensitive to the chain stiffness. 

In order to quantify the amount of adsorbed monomers as a function of the distance z 
from the charged wall, we have also studied N*(z) as defined by Eq. (jHJ). Our results are 
reported in Fig. EJ At A e = [see Fig. 11(a)] the fraction of adsorbed monomers N*(z) is 
always larger with growing N m even near the interface. On the other hand, at A e = 0.951 
[see Fig. 0(b)] : (i) N*(z) gets smaller with growing N m near the interface (roughly for 
z/a < 0.65) and (ii) N*(z) is considerably reduced compared to the A e = 0-case. For 
instance (with N m = 12) at z/a = 0.9 (corresponding to a layer thickness at A e = 0.951), 
N*(z) can be as large as 0.4 for A e = against only 0.2 for A e = 0.951. 

2. PE orientation 

To gain further insight into the properties of rod-like PE adsorption, we have plotted S(z) 
as given by Eq. ljllj) in Fig. El so as to characterize the PE orientation with respect to the 
charged interface. At A e = 0, Fig. E(a) shows that in the vicinity of the interface (roughly 
for z < a) that the rod-like PEs tend to be parallel to the interface-plane with growing N m , 
i.e., S(z) — > — 1/2 (see also Fig. EJ. This effect is obviously due to the electrostatic wall-PE 
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FIG. 5: Profiles of S(z) for different chain length N m with oqL 2 = 64 (systems A — E). (a) A e = 0. 
(b) A t = 0.951. 

binding whose strength increases linearly with N m . For the longest chains (N m = 12) non- 
negligible positive S'(z)-values are reported at intermediate distance from the wall (roughly 
for 2 < z/a < 6) signaling a non-trivial orientation correlation with respect to the interface- 
plane that will be properly discussed later. Sufficiently away from the wall, the rod-like PEs 
are randomly oriented [i. e., S(z) — > 0] as it should be. 

In the presence of image forces [A e = 0.951 - see Fig. Efb)] the S(z) behavior is more 
complex. A comparison with Fig. Efa) corresponding to A e = immediately shows that 
repulsive image forces tend to inhibit the alignment of the rod-like PEs with respect to 
the interface-plane near contact. This effect will be especially vivid at lower surface charge 
density Co, as we are going to show later. 

The non-monotonic behavior of S(z) near contact at A e = 0.951, similar to that reported 
for n(z) in Fig. ^b), is the result of two antagonistic forces: (i) the repulsive image driving 
force that scales like 1/z and (ii) the attractive wall-monomer one that scales like z. As in 
the case with A e = 0, (i) the degree of PE-wall parallelism increases with growing N m near 
contact and (ii) far enough from the wall the PEs are randomly oriented. 



3. Fluid charge 

Another interesting property is provided by the net fluid charge parameter o~*(z) [Eq. (|13j)] 
that describes the screening of the charged interface. The profiles of o~*(z) for different N m 
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FIG. 6: Profiles of the reduced net fluid charge c*{z) for different chain length N m with a^L 2 = 64 
(systems A - E). (a) A e = 0. (b) A e = 0.951. 

can be found in Fig. |3 At A e = [see Fig. EJa)], it is shown that for long enough chains 
(here N m > 4) the substrate gets locally overcharged as signaled by o~*(z) > 1. Physically, 
this means that the (integrated) local charge of the adsorbed monomers^ is larger in absolute 
value than that of the substrate's surface charge. In other words, the plate is overscreened by 
the adsorbed PE chains. Fig. Elf a) indicates that the degree of overcharging increases with 
N m as expected from the behavior of multivalent counterions. Upon inducing polarization 
charges [A e = 0.951 - see Fig. Efb)] overscreening is maintained and weakly disturbed, 
proving that the latter is robust against repulsive image forces. 

B. Influence of substrate's surface-charge density 

To complete our investigation, we would like to address the influence of the substrate 
charge density on the PE adsorption. In this respect, we consider (at fixed N m = 8) two 
additional values for the surface charge density: <JqL 2 = 32 and 128 corresponding to the 
systems F and G, respectively (see Table EJ). 

1. Monomer distribution 

The plots of the monomer density n(z) for various values of <JqL 2 can be found in Fig. 
[7| Typical microstructures of systems F and G are sketched in Fig. |H1 At A t = [see Fig. 
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FIG. 7: Profiles of the monomer density n{z) for different parameters of surface-charge density 
oqL 2 with N m = 8 (systems D, F,G). The case o~qL 2 = 64 (system D) from Fig. reported here 
again for easier comparison, (a) A e = 0. (b) A e = 0.951. 

Ufa)], the monomer density at contact increases with cr as it should be. Interestingly, the 
local maximum in n(z), present at small cq (here cr L 2 < 64), vanishes at large <7o [see Fig. 
Ufa)]. This feature is the result of a ao-enhanced driving force of adsorption that overcomes 
chain-entropy effects at large enough oq. The strong adsorption at <JqL 2 = 128 leads to a 
flat PE layer as well illustrated in Fig. |S| 

By polarizing the interface (A e = 0.951), it can be seen from Fig. Ulb) and the snapshot 
from Fig. [SJthat there is a strong monomer depletion near the interface for a^L 2 = 32. This 
feature is due to the combined effects of (i) image-PE repulsion and (ii) a lower electrostatic 
wall-PE attraction. Upon increasing <7o the monomer density near contact becomes larger, 
and concomitantly, the location of the maximum in n(z) is systematically shifted to smaller 
z. It is to say that the thickness of the adsorbed PE layer decreases with cq. We also 
expect that, at very large <To (not reported here), this maximum vanishes leading to a purely 
attractive effective wall-PE interaction. 



2. PE orientation 



At A e = 0, Fig. Efa) shows that near the charged interface (when z < a) the degree of 
parallelism between the rod-like PE and the interface increases with growing ctq as indicated 
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FIG. 9: Same as Fig. [3 but for S(z). 

by S(z) — > —1/2. This observation is merely due to the electrostatic wall-PE binding whose 
strength scales like <7o at fixed N m . 

In the presence of image forces (A e = 0.951), Fig. Efb) demonstrates again for a^L 2 = 
32 [see also Fig. Efb) for comparison)] a strongly non-monotonic behavior of S(z) near 
the interface. This feature is fully consistent with the ideas that (i) image forces become 
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FIG. 10: Same as Fig. Qbut for a*(z). 

especially strong at low o"o and (ii) repulsive image forces induce orientational disorder as 
previously established. This finding leads to the important general statement: 

• Repulsive image forces at low surface charge density induce orientational disorder near 
the interface. 



3. Fluid charge 

The profiles of o~*(z) for different o§L 2 can be found in Fig. EH At A e = [see Fig. 
[TUT a)], it is found that the planar interface gets always locally overcharged as signaled by 
o~*(z) > 1. The location of the maximum in cr*(z) is shifted to lower z with increasing o~q. 

Upon inducing polarization charges [A e = 0.951 - see Fig. ITOT b)] overscreening is still 
there. However, at o~qL 2 = 64, there is a non-negligible shift of the maximum of about 0.5a. 
The distance at which the substrate is compensated [i.e., where o~*(z) = 1] corresponds to 
1.72a (2.54a) for A e = (A e = 0.951) leading to a neutralization z-shift of 0.72a. 

IV. SUMMARY 

To conclude, we have performed MC simulations to address the behavior of rod-like poly- 
electrolytes at oppositely charged planar surfaces. The influence of image forces stemming 
from the dielectric discontinuity at the charged interface was also analyzed. We have con- 
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sidered a finite and fixed monomer concentration in the dilute regime for relatively short 
chains. 

In the absence of image forces (A £ = 0), our main findings can be summarized as follows: 

• At moderately charged interfaces, only (very) short rod-like PEs (here 2 < N m < 8) 
experience a short-ranged repulsion near the interface. For longer rod-like PEs the 
effective wall-PE interaction becomes purely attractive. This behavior is in contrast to 
that occurring with flexible PEs, where the chain-entropy is larger leading to stronger 
entropy-driven depletion. 

• Near the charged interface, the rod-like PEs get more and more parallel to the interface- 
plane when the chain length N m is increased. Concomitantly, the strength of the 
substrate-charge reversal is A m -enhanced. 

• Upon increasing the substrate-surface-charge density o"o it was demonstrated that: (i) 
The monomer adsorption is enlarged and the wall-PE effective interaction becomes 
purely repulsive for high enough cr - (h) The degree of parallelism (near the interface) 
between the interface-plane and the rod-like PE is enhanced. 

The main effects stemming from repulsive image forces (A e = 0.951) as revealed by this 
study follows: 

• The monomer adsorption is reduced at finite A e and the n(z)-profiles become similar 
to those obtained with flexible chains, in contrast to what was reported at A e = 0. 

• Repulsive image forces induce PE orientational disorder near the interface. 

• The substrate-charge reversal is robust against repulsive image forces. 
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